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w ' In Ref. [l] we reported an estimate of the critical exponent for the divergence of the 

localization length at the quantum Hall transition that is significantly larger than those 
reported in the previous published work of other authors. In this paper, we update our 
finite size scaling analysis of the Chalker-Coddington model and suggest the origin of 

►^ ' the previous underestimate by other authors. We also compare our results with the 

— ^ , predictions of Liitken and Ross.l^ 
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1. Introduction 



When a strong magnetic field is appfied perpendicular to an ideal two dimensional 
electron gas the kinetic energy of the electrons is quantized according to the formula 
k> ' En = {n-\-l/2)Huj. Here, n is a non-negative integer and to is the cyclotron frequency 

^ . u! = eB/ra. These Landau levels are highly degenerate and the density of states 

becomes a series of equally spaced delta functions. This degeneracy is broken by 
disorder and the Landau levels are broadened into Landau bands. Most of the 
electron states are Anderson localized with the exception of the states at the center 
of the Landau level where the localization length ^ has a power law divergence 
described by a critical exponent v 

i^\E-E,\-'' . (1) 

When the Fermi level is in a region of localized states, the Hall conductance is 
quantized in integer multiples of e^/ft. This effect is known as the quantum Hall 
effect. Transitions between consecutive quantized values occur when the Fermi 
level passes through the center of a Landau band. This is a quantum phase tran- 
sition. It is characterized by a two critical exponents. One is the critical exponent 
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V mentioned above and the other is the dynamic exponent z, which describes tem- 
perature dependence. 

The quantum Hall transition has been the subject of careful experimental study. 
The inverse of the product of the critical and dynamic exponents, k = 1/vz, has been 
measured very precisely (see Table [1]). The value of this product appears to be quite 
universal; the same value has been obtained in measurements in an Ala;Gai_a;As 
heterostructure^ and in graphene.l^The main problem is that an independent mea- 
surement of the dynamic exponent is needed to disentangle the values of the two 
exponents and, unfortunately, this has been measured much less precisely. 

The quantum Hall transition has also been the subject of numerous numerical 
studies in models of non-interacting electrons. In earlier work, a consensus was 
reached (see Table[2]) that v « 2.4 in apparent agreement with experiment. However, 
in 2009, we published^J a numerical analysis of the Chalker-Coddington model in 
which we found a value of the exponent that was about 10% larger; a result which 
has since been confirmed by other authors (see Table [3]). There now seems to be a 
consensus that the previous numerical work underestimated the exponent and the 
apparent agreement with experiment was a coincidence of errors. Below we discuss 
the reason for the previous underestimate of the exponent. 

Table 1 . Experimental values of critical exponents for the quantum Hall 
transition. 

K = 1/ ZV V Z 

Experiment of Li et al.E! (Al^jGai-:! As) 0.42 ± .01 fs 2.38 « 1 

Experiment of Giesbers et al.l^ (graphene) 0.41 ± .04 



Table 2. Earlier estimates of the critical exponent v. 



Chalker and CoddingtorF-J 


2. 5 ±.5 


Huckestein and Kramejl^l 


2.34 ± 


.04 


MieclP 


2. 3 ±.08 


HuckesteiiilQJ 


2.35 ± 


.03 


Huo and BhatlPJJ 


2.4±.l 


Lee and Wan^ 


2.33 ± 


.03 


Cain et al.1121 


2. 37 ±.02 









Table 3. Recent estimates of the critical exponent i/. 

Slevin and Ohtsukff 2.593(2.587, 2.598] Obuse et al.lJ^ 2.55 ± .01 
Dahlhaus et alllSI 2.576 ± .03 Amado et al.E! 2.616 ± .014 



Another issue concerns the value of the dynamic exponent. For models of non- 
interacting electrons the dynamic exponent is known exactly, z — 2. Howe ver , this 
value cannot be compared directly with the experiment. Burmistrov et al.'^'^ have 
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emphasized the distinction between the different dynamical exponents that occur 
in the problem. In the experiment of Li et alP it seems clear that the dynamic 
exponent that is being measured describes the divergence of the phase coherence 
length on approaching zero temperature 

£^ ~ T-^ . (2) 

It also seems reasonably safe to suppose that electron-electron interactions are the 
source of the electron dephasing. This does not mean however that electron-electron 
interaction are relevant in the renormalization group (RG) sense and that the quan- 
tum Hall transition is described by a fixed point in a theory of interacting electrons. 
A clear discussion of this can be found in Ref. [T71 Also, as described in Ref. [T51 
it is thought that short range interactions are irrelevant in the RG sense and that 
only long range Coulomb interaction are relevant and would drive the system to a 
different interacting fixed point. 

It has been pointed out to us by Alexei Tsvelik, that the value of the exponent 
we have found for the Chalker-Coddington model is very close to that predicted by 
Liitken and Ross. In a series of papers (Ref. |2]and references therein) these authors 
have argued that modular symmetry strongly constrains the possible critical theories 
of the quantum Hall transition. While we cannot claim to understand the details 
of the theory of Liitken and Ross, we attempt below to compare some of their key 
predictions with the results of our finite size scaling analysis. 

2. Method 

We calculated the Lyapunov exponents of the product of the transfer matrices for 
the Chalker-Coddington model.^This model describes electron localization in a two 
dimensional electron gas subject to a very strong perpendicular magnetic field. The 
basic assumption is that the random potential is smooth on the scale of the magnetic 
length. The electron wavefunctions are then concentrated on equipotentials of the 
random potential with tunneling between equipotentials at saddle points of the 
potential. In the Chalker-Coddington model this system is modeled by a network 
of nodes and links. A parameter x, which is essentially the energy of the electrons 
measured in units of the Landau band width relative to the center of the Landau 
band, fixes the tunneling probability at the nodes. A random phase distributed 
uniformly on [0, 2tt) is attached to each link to refiect the random length of the 
contours of the potential. For further details we refer the reader to the original 
article of Chalker and CoddingtorM-' and to the more recent review by Kramer et 
allEI 

We considered the transfer matrix product associated with a quasi-one dimen- 
sional geometry with N nodes in the transverse direction and L nodes in the lon- 
gitudinal directionO The Lyapunov exponents of this random matrix product were 

''For the detailed formulae see Ref. [l] 
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Fig. 1. The product T of the smallest positive Lyapnuov exponent 7 and the number of 
nodes in the transverse direction Af = 4, 6, 8, 12, 16, 24, 32, 48, 64, 96, 128, 192, 256 for the Chalker- 
Coddington model. The error in the data is much smaller than the symbol size. The lines are even 
order polynomial fits to the data for each N . 



estimated using the standard method! I I The Lyapunov exponents are defined 
by taking the Umit L -^ 00. By truncating the matrix product at a finite L, an esti- 
mate of the Lyapunov exponents was obtained. The sample to sample fluctuations 
of this estimate decrease with the inverse of the square root of L. We performed 
a single simulation for each pair of x and N and truncated the transfer matrix 
product at a value of L that allowed estimation of the smallest positive Lyapunov 
exponent 7 with a precision of 0.03%, except for the largest values of A^ = 192 
and 256 where the precision was relaxed to either 0.05% or 0.1%. To ensure that 
simulations for different pairs of x and N were independent, we used the Mersenne 
Twister pseudo-random number generator MT2203 of Matsumoto et al.!^ provided 
in the Intel Math Kernel Library. All the simulations used a common seed. Indepen- 
dence was ensured by the use of a unique stream number for each simulation. We 
imposed periodic boundary conditions in the transverse direction for which choice 
the Lyapunov exponents are even functions of x. It is known that there is a critical 
point at the center of the Landau band, x = 0, and that, when the Fermi energy is 
driven through this point, the transition between Hall plateaux occurs. 

To extract estimates of critical exponent and other quantities we used finite size 
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Fig. 2. The same data as in Fig. [T] but focussing on the critical point at x 
variation of F with A^ at a; = is due to irrelevant scaling variables. 

scalingcl In this method, the behavior of the dimensionless quantity 

T{x,N) =-fN , 



0. The residual 



(3) 



is analyzed as a function of both x and TV. In the absence of any corrections to 
scahng we would expect this behavior to be described by the following finite size 
scaling law 



where i^ is an a priori unknown but universal scaling function and 

a — l/i/ . 



(4) 



(5) 



Note that we have imposed the condition that F must be an even function of x. 

The actual behavior of F as a function of x for different N is shown in Fig.[T]and, 
in more detail around a; = 0, in Fig. [2l (The lines in the figures are polynomial fits. 
They will be discussed below.) According to Eq. ^ curves for different N should 
have a common crossing point at x — 0. However, it is clear from Fig. [5] that this is 
only approximately correct and that F is not exactly independent of iV at a; = but 
varies by a several percent over the range of N studied. These corrections to scaling 



^This metho d was first applied to Anderso n locali zation at about the same time by Pichard and 
Sarm a I I and by MacKinnon and Krame r I I See Ref. 1261 for a pedagogical discussion. 



February 25, 2013 10:53 WSPC/INSTRUCTION FILE 

31evin-Ohtsuki'Localisation'2011 



6 KEITH S LEVIN, TO MI OHTSUKI 

arise because of the presence of irrelevant scaling variables. These are variables with 
negative scaling exponents. Their effect is negligible for large N but their presence 
may lead to significant corrections at small N . This is consistent with what we see 
in Fig. H 

To take account of corrections to scaling, we follow Ref. 27 and Ref. 28 and 
generalize the finite size scaling law 

T ^ F {N^''vQ{x),Ny'vi{x),Ny-V2{x),- ■ ■) . (6) 

In Eq. ([6]) , Wo is the relevant scaling variable and wi , W2 , ■ • ■ are the irrelevant scaling 
variables. The associated exponents j/i, 2/2i ■ ' ' ^-^e negative. The inclusion of irrele- 
vant corrections permits the residual N dependence at x = seen in Fig. [2] to be 
modeled. This form also allows for additional corrections due to non-linearities of 
the scaling variables as functions of a;. To impose the condition that F must be an 
even function of x we restrict all the scaling variables to be even functions of x. In 
addition, since the critical point is at a; = 0, we impose the condition that wo(0) = 0. 
To fit the data, the function F is expanded as a Taylor series in all its arguments, 
and similarly the scaling variables. The coefficients in the Taylor series, together 
with the various exponents, play the role of fitting parametersU The orders of trun- 
cation of the Taylor series are chosen sufficiently large to obtain an acceptable fit 
of the data (as measured using the x^-statistic and the goodness of fit probability). 
We have attempted this procedure with both one and two irrelevant corrections. 
Unfortunately, a stable fit of the data has eluded us. We have found that several 
fits of the data are possible. However, these fits do not yield mutually consistent 
estimates of the critical exponent. 

3. Rudimentary finite size scaling 

To circumvent the difficulties described in the previous section we resorted to a 
less sophisticated approach in which we abandoned the attempt to fit all the data 
in a single step. Instead, we fitted the data for each N independently to an even 
polynomial of x. For each N the order of the polynomial was chosen just large 
enough to give an acceptable goodness of fit. For TV = 4, a quadratic was sufficient, 
while for N = 256 a sixth order polynomial was required. From these polynomials 
we estimated the curvature C of F at a; = 0. The results are plotted in Fig. [3] 
The precision of the estimation of the curvature varies between 0.07% and 0.28%. 
According to (j4|) the curvature at x = should vary with iV as a power law, 



dx^ 



a:=0 



This, of course, neglects corrections to scaling due to irrelevant variables and, indeed, 
a straight line fit to all the data does not yield an acceptable goodness of fit. 

'^Some extra conditions on tiic coefEcients must be imposed to ensure that the model to be fitted 
to the data is unambiguous. 
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N 



Fig. 3. The curvature C at x = (see Eq. iQ), obtained from the polynomial fits shown in Figure 
[T] plotted as a function of A'^. The solid line is a straight line fit to the data for the largest five 
values of N. The slope corresponds to u = 2.607. The dashed line is a straight line with slope 
corresponding to i^ = 2.34 and passing though the TV = 4 data point. The error in the data for the 
curvature is much smaller than the symbol size. 

However, if data for iV < 48 are excluded, acceptable goodness of fits are obtained. 
The estimates of the critical exponent obtained in this way are tabulated in Tabled 
In Fig. [3] we have also plotted two lines. One is a solid line that corresponds 
to the straight line fit for 64 < A^ < 256. The second is a dashed line with a 
slope corresponding to the estimate of HuckesteirR-^ of the critical exponent and 
passing through the datum for the curvature at A^ = 4. While the main effect of the 
irrelevant corrections is the iV dependent shift in the ordinate that is clearly visible 
in Figl21 a smaller but not negligible effect on the curvature is also apparent in Fig|31 
In our opinion, this is the reason why the critical exponent was underestimated in 
previous work. The precision of the numerical data was insufficient, and the range 
of N considered too small, for the irrelevant correction to be properly taken into 
account. 



4. Comparison v^rith the predictions of Liitken and Ross 

The predictions of the theory of Liitken and Ross that can be compared with the 
present work are for the critical exponent and the leading irrelevant exponent. Their 
theory contains a single unknown parameter, the central charge c. In terms of this 
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Table 4. Estimates of the critical exponent v obtained from 
linear fits of InC versus \nN . 

range of A'' considered v Confidence intervals (±20-) 

64 < Af < 256 2.607 [2.598, 2.615] 

96 < Af < 256 2.599 [2.583, 2.616] 

128 < Af < 256 2.590 [2.586,2.614] 

parameter they predict that the critical exponent iu 

10.42050633345819 • • ■ 
.= ^ ■ (8) 

In addition they predict that the leading irrelevant exponent and the critical expo- 
nent are related by 

y = -a= -1/v . (9) 

While we are not aware of any justification for this, assuming c — A gives 

i/ = 2.60512- •• (10) 

The most precise estimate in Table |4] is 

1/ = 2.607 ± 0.004 . (11) 

In fact, all of the estimates in Table |4] are consistent with the Liitken and Ross 
value. Turning to the irrelevant exponent, the situation is, unfortunately, much less 
clear. In Table [5] we tabulate some previous estimates of the irrelevant exponent. 
The Liitken and Ross prediction is 

2/ = -0.383859 •• • (12) 

The estimate of HuckesteirP^ is consistent with this but subsequent estimates by 
Wang et al.l^^ are somewhat ambiguous. In our opinion, further confirmation is 
needed before reaching a conclusion. 

Table 5. Estimates from the literature of the leading irrelevant exponent y and the 
quantities used to estimate it. 



Huckestein[2ZI 


y = 


-0.38 ± .04 


Lyapunov exponents 


Wang et all221 


y~ 


-0.52 


geometric average of the two-terminal conductance 


Wang et al.ESI 


y~ 


-0.72 


arithmetic average of the two-terminal conductance 



As mentioned above an unambiguous fit using Eq. ([6]) has not proved possible. 
The best we can do at present is to check the consistency of the Liitken and Ross 
values with our data by fixing both v and y to these values when fitting. A series 

"^We noticed an error (a factor of 2) in Ref. (2] We thank Graham Ross for confirming this. 
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Table 6. Fits to the 379 data points shown in Figure [T] 
rii number of parameters x^ goodness of fit Vc 



3 3 16 398 0.10 0.807 ± .0005 

3 4 20 347 0.66 0.804 ± .0015 

3 5 24 343 0.67 0.801 ± .0035 



Table 7. Estimates of Fc obtained using published 
estimates of oo and Eq. I I15II 

ao F, 

Obuse et al.ESl 2.2617 ± 0.0006 0.8222 ± .0019 
Evers et al.EU 2.2596 ± 0.0004 0.8156 ± .0013 



of such fits are tabulated in Table [51 In the fits only a single irrelevant variable is 
assumed and non-linearities in the scaling variables are ignored, i.e. 

VO = Vq2X^ , vi = Wio . (13) 

The scaling function is expanded as a Taylor series to order tiq in the relevant field 
and order ni in the irrelevant field. One of the important quantities that can be 
estimated in this way is Fc, which is defined by 

Fc= lim r(2; = 0,7V) = i^ (0,0,...) . (14) 

The quantity is significant because, if the quantum Hall critical theory has conformal 
symmetry, it is related to the multi-fractal exponent ao that occurs in the multi- 
fractal analysis of the wavefunction distribution at the critical point by 

F, = 7r(ao-2) . (15) 

Some estimates of Fc obtained using published estimates of a^ and Eq. (|15p are 
tabulated in Table [71 Our numerical estimate of Fc is not completely consistent 
with those in the table. The reason for this is not yet clear. 

5. Discussion 

The agreement between our estimate for the critical exponent and prediction of 
Liitken and Ross is tantalizing but is it accidental, or does it have a deeper sig- 
nificance? A more precise numerical estimate of the irrelevant exponent is clearly 
highly desirable. In addition, the prediction of Liitken and Ross for the flow diagram 
in the {(Jxy,<7xx) plane should also be amenable to numerical verification. We also 
need to know if a central charge c = 4 is physically justified. 

Quite apart from whether or not the Liitken and Ross theory is exact for non- 
interacting electrons, the important question remains of clarifying the role of the 
electron-electron interactions in the observed critical behavior at the quantum Hall 
transition. More work along the lines of Ref. 1321 might be very helpful in this regard. 



February 25, 2013 10:53 WSPC/INSTRUCTION FILE 

31evin-Ohtsuki'Localisation'2011 



10 KEITH SLEVIN, TOMI OHTSUKI 

Acknowledgments 

We would like to thank Alexei Tsvelik for bringing the work of Liitkcn and Ross 
to our attention. This work was supported by Grant-in- Aid 23540376 and Korean 
WCU program Project No. R31-2008-000-10059-0. 

References 

1. K. Slevin and T. Ohtsuki, Physical Review B 80, 041304(R) (2009). 

2. C. A. Lutken and G. G. Ross, Pliysics Letters B 653, 363 (2007). 

3. K. V. Klitzing, G. Dorda, and M. Pepper, Pliysical Review Letters 45, 494 (1980). 

4. D. Yoshioka, The quantum Hall effect, Springer series in solid-state sciences, 
(Springer, Berlin ; New York, 2002). 

5. W. Li et al, Physical Review Letters 102, 216801 (2009). 

6. A. J. M. Giesbers et al., Physical Review B (Condensed Matter and Materials Physics) 
80, 241411 (2009). 

7. J. T. Chalker and P. D. Coddington, Journal of Physics C: Solid State Physics 21, 
2665 (1988). 

8. B. Huckestein and B. Kramer, Physical Review Letters 64, 1437 (1990). 
. Mieck, EPL (Europhysics Letters) 13, 453 (1990). 
. Huckestein, EPL (Europhysics Letters) 20, 451 (1992). 
. Huo and R. N. Bhatt, Physical Review Letters 68, 1375 (1992). 
.-H. Lee and Z. Wang, Philosophical Magazine Letters 73, 145 (1996). 

Cain, R. A. Romer, and M. E. Raikh, Physical Review B 67, 075307 (2003). 
. Obuse et al.. Physical Review B 82, 035309 (2010). 

P. Dahlhaus, J. M. Edge, J. Tworzydlo, and C. W. J. Beenakker, Physical Review 
B 84, 115133 (2011). 

Amado et al. Physical Review Letters 107, 066402 (2011). 
S. Burmistrov et al., Annals of Physics 326, 1457 (2011). 
M. M. Pruisken, International Journal of Modern Physics B 24, 1895 (2010). 
Kramer, T. Ohtsuki, and S. Kettemann, Physics Reports 417, 211 (2005). 
Shimada and T. Nagashima, Progress of Theoretical Physics 61, 1605 (1979). 
MacKinnon and B. Kramer, Zeitschrift fiir Physik B Condensed Matter 53, 1 
(1983). 
22. M. Matsumoto and T. Nishimura, in Monte Carlo and Quasi-Monte Carlo Methods 
1998, edited by H. Niederreiter and J. Spanier (Springer, Berlin ; New York, 1999), 
pp. 56-69. 

L. Pichard and G. Sarma, Journal of Physics C: Solid State Physics L127 (1981). 
L. Pichard and G. Sarma, Journal of Physics C: Solid State Physics L617 (1981). 
. MacKinnon and B. Kramer, Physical Review Letters 47, 1546 (1981). 
. J. Amit and V. Martin-Mayor, Field theory, the renormalization group, and critical 
phenomena, 3rd ed. ed. (World Scientific, New Jersey ; London, 2005). 
Huckestein, Physical Review Letters 72, 1080 (1994). 
Slevin and T. Ohtsuki, Physical Review Letters 82, 382 (1999). 
Wang, Q. Li, and C. M. Soukoulis, Physical Review B 58, 3576 (1998). 
Obuse et al. Physical Review Letters 101, 116802 (2008). 

Evers, A. Mildenberger, and A. D. Mirlin, Physical Review Letters 101, 116803 
(2008). 
32. B. Huckestein and M. Backhaus, Physical Review Letters 82, 5100 (1999). 



9. 


B 


10. 


B 


11. 


Y 


12. 


D 


13. 


P. 


14. 


H 


15. 


J. 




B 


16. 


M 


17. 


I. 


18. 


A 


19. 


B 


20. 


I. 


21. 


A 



23. 


J. 


24. 


J. 


25. 


A 


26. 


D 


27. 


B 


28. 


K 


29. 


X 


30. 


H 


31. 


F. 



